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PREFACE 

This  paper  is  based  on  research  performed  for  the  United 
States  Air  Force  under  Project  RAND  and  has  been  reported  in 
other  forms.  Any  views  expressed  in  this  paper  should  not  be 
interpreted  as  reflecting  the  official  opinion  or  policy  of  any  of 
The  RAND  Corporation's  governmental  or  private  research 
sponsors. 

This  paper  was  prepared  for  publication  in  11  Nuovo  Cimento. 


ON  THE  FOUNDATIONS  OF  RELATIVISTIC  ENERGY  MECHANICS 

by 

Dominic  G.  B.  Edelen 
The  RAND  Corporation 
Santa  Monica,  California 

SUMMARY 

The  Einstein  theory  of  general  relativity  is  shown  to  yield  a  general  me¬ 
chanics  of  continuous  media  under  the  assumption  that  the  momentum-energy 
tensor  admits  a  unique  time-like  eigenvector.  Physical  interpretations  of  the 
governing  equations  are  derived,  together  with  constitutive  relations  for  gener¬ 
al  and  isotropic  materials.  It  turns  out  that  the  mechanics  can  always  be  viewed 
as  describing  the  flow  of  rest-energy.  Invariant  requirements  for  the  existence 
of  a  stress  potential  are  obtained,  the  satisfaction  of  which  leads  to  a  decom¬ 
position  and  partial  evaluation  of  the  rest-energy.  The  Einstein  field  equations 
are  shown  to  imply  the  existence  and  uniqueness  of  an  intrinsic  energy  density 
for  any  material  medium  (intrinsic  Immutable  mass).  The  usual  procedure  of 
adding  conditions  to  the  Einstein  theory  in  order  to  obtain  an  analogous  intrinsic 
quantity  is  thus  unnecessary.  The  path  density,  which  defines  the  intrinsic 
energy,  is  shown  to  be  path  independent  in  an  appropriate  sense  if  the  stresses 
admit  a  stress  potential.  This  suggests  a  decomposition  of  the  generalised 
stresses  and  leads  to  a  fundamental  differential  relation  on  the  trajectories  of 
the  energy  flux.  Natural  definitions  of  intrinsic  temperature  and  intrinsic  en¬ 
tropy  density  are  direct  consequences  of  the  fundamental  differential  relation 
and  lead  to  generalized  thermodynamic  descriptions  which  include  the  effects 
of  gravitational  radiation. 


1.  INTRODUCTION 


Sine*  the  promulgation  of  the  general  theory  of  relativity,  science  has  come 
to  view  the  problem  of  describing  general  media  as  equivalent  to  the  problem  of 
stating  the  momentum-energy  tensor  for  such  media.  For  certain  special  cases, 
the  appropriate  forms  for  the  momentum-energy  tensor  are  known.  These  have 
been  most  useful  in  gleaning  the  dynamic  interplay  of  geometry  and  physics. 
Unfortunately,  without  significant  extension,  these  cases  regulate  the  investi¬ 
gator  to  the  consideration  of  only  the  simplest  of  possible  physical  models.  In 
view  of  the  fundamental  nature  of  the  Einstein  theory,  it  should  be  possible  to 
develop  a  simple  and  consistent  mechanics  of  general  systems  which  involves  a 
minimum  of  assumptions  and  which  yields  direct  physical  interpretations  of  all 
quantities  involved.  The  purpose  of  this  paper  is  to  lay  the  requisite  foundations 
for  such  a  mechanics. 

2.  THE  MOMENTUM-ENERGY  TENSOR  OF  A  MATERIAL  MEDIUM 

Let  ^  be  an  Einstein— Riemann  space*  with  metric  tensor  h^Q  and 

associated  momentum-energy  tensor  T Capital  Latin  indices  are  assumed 

to  range  through  the  integers  from  one  to  four.  A  semicolon  will  be  used  to 

denote  covariant  differentiation,  the  comma  being  reserved  to  denote  partial 

differentiation.  The  usual  notation  for  mixing  and  alternating  will  be  used. 

Following  the  line  of  thought  introduced  by  Synge**,  it  appears  natural, 

and  indeed  necessary,  to  require  that  any  region  of  which  is  filled  by  a 
- • - 

A  four-dimensional  metric  space  of  the  hyperbolic-normal  type  in  which 
Einat^g's  field  equations  are  satisfied.  The  signature  of  £  is  therefore  —2. 

Synge,  J.  L. ,  Relativity:  The  General  Theory,  North-Holland  Publishing 
Co.,  Amsterdam  (1960),  pp.  leg— 169.  "" 


material  medium  be  such  that  its  associated  momentum-energy  tensor  admit  a 
unique  time-like  eigenvector.  The  precise  characterization  is  as  follows: 
DEFINITION  1.  A  symmetric  tensor  field  T^  is  a  material  momentum-energy 
tensor  associated  with  a  region  ^  of  ^  if  and  only  if 

(2.1)  TAB  n  -  0 

,o 

A 

and  there  exists  a  unique,  time-tike,  unit  vector  field  W  and  a  nonzero  scalar 
p  on  Z5  such  that 

(2.2)  TAWA-pWB. 

Careful  attention  should  be  paid  to  the  requirement  that  the  vector  field 

A  » 

W  be  unique.  This  uniqueness  states  that  any  point  P  of  the  region  &  has 

A 

associated  with  it  a  unique  time-direction  defined  by  W  ,  end  hence  the  material 
medium  filling  has  a  uniquely  defined  time-orientation.  Without  such  a  result, 
it  is  difficult,  if  not  altogether  impossible— to  speak  meaningfully  of  the  simplest 
everyday  properties  of  what  we  are  accustomed  to  think  of  as  material  bodies. 

Care  must  also  be  exercised  in  the  use  of  Def.  1,  for  we  have  to  distinguish  be¬ 
tween  the  various  cases  which  may  arise. 

DEFINITION  2.  A  material  momentum-energy  tensor  is  said  to  be  of  class 

m(m  •  1,  2,  3,  4)  if  and  only  if  the  multiplicity  of  the  eigenvalue  n  Is  m. 

\ 

In  order  to  guarantee  the  unlclty  of  the  eigenvector  W  ,  material  momentum- 

energy  tensors  of  class  m  must  be  such  that  their  characteristics  (referred  to 

Jordan  normal  form)  commence  with  the  integer  m.#  Since  space-Uke  and 

*For  a  detailed  discussion  of  such  characteristics  see  Hlavat/,  V. , 
"Contribution  to  the  Theory  of  General  Geometrodynamics"  (to  be  polished). 


tlme-like  quantities  usually  differ  markedly  in  their  numerical  values  when 
calculated  with  other  than  what  may  be  termed  exotic  units,  it  is  natural  to 
expect  that  material  momentum-energy  tensors  of  class  one  predominate. 
For  this  reason,  we  confine  our  attention  in  the  remainder  of  this  paper  to 
such  momentum-energy  tensors. 


3.  CHARACTERIZATION  OF  MATERIAL  MOMENTUM-ENERGY  TENSORS 
OF  CLASS  ONE 


The  definitions  given  in  the  previous  section  are  implicit,  for  they  do 

AB 

not  exhibit  the  detailed  form  and  nature  of  T  .  This  section  is  devoted  to 
developing  the  explicit  characterization  which  is  needed  in  order  to  obtain 
physical  interpretations. 

Consider  the  idempotent  tensor  field 

<3*i>  ■5--;-wawb 

on  .  With  it  we  may  define  a  projection  operator  w  over  the  set  of  all 

geometric  objects  on  zV  as  follows: 

,*A. » . .  A  D  ,C ... 

,^IB.  .  .  9C’  ‘  **B* '  ‘u. . . 


In  particular,  we  have 

«•»  "*A;B,-W(A;B)-W<a'V+  W(A.B]  *W[A*B]' 

where 

(3.3)  ^Ad,fwA;BwB- 

(Here  we  have  used  the  fact  that  WA  ■  1  implies  W A  wb.q  ■  0> )  A  direct 
application  of  well  known  techniques  together  with  definitions  1  and  2  then  yields 
the  following  result.  The  most  general  symmetric  tensor  field  of  class  one 
which  satisfies  (2. 2)  is  given  by 


(3.4) 


TAB.|lWAWB+CAB, 


AB 


where  a  is  any  symmetric  tensor  field  such  that 

(3.5) 
end 


.  AB.  AB 
»(o  )  «  a 


(3.6)  aABVB-^VA 
if  end  only  if  V&  -  0. 

With  the  above  result,  all  that  remains  in  order  for  T  to  be  material 
is  satisfaction  of  equations  (2. 1).  Substituting  (3.4)  Into  (2. 1)  gives 

(3.7) 


,AB 


P  WA  +  WA  (*WB).b  ♦  0ABb  -  0. 


A  AB 

Now,  by  assumption  we  have  W  »  1,  while  (3.5)  implies  that  <7  W} 

It  thus  follows  that  WA  ■  0  and 

<jAB;B  WA  “  "aAB  WA;B  "  “ffAB  W(A;B) 

AB  AB 

since  a  is  symmetric.  Noting  that  o  W(A  B)  a  scalar  and  that  o 


0. 


AB 


satisfies  (3.5),  we  have 
AB 


AB 


°  *  AB' 


(3.6) 

where  (vid  3.2) 

(3.9) 

is  the  Born  rate-of-etrain  tensor  associated  with  the  vector  field  W “ .  * 
Hence,  if  we  contract  (3.  7)  with  we  are  led  to  the  scalar  equation 

(3.10) 


.  ,„B.  AB 

}-*m°  *AB’ 


;B 

With  the  aid  of  the  Lie  derivative,**  we  can  write  (3. 10)  in  the  more  suggestive 


Synge;  op.  clt. .  p.  172,  eq.  62.  The  difference  in  sign  between  our 
statement  and  that  of  Synge  arises  because  Synge  assumes  a  signature  +2  while 
we  have  assumed  a  signature  -2.  The  physical  interpretation  of  « .  _  is  die- 
cuseet^n  Sec.  7. 

Yano.  K.,  The  Theory  of  Lie  Derivatives  and  Its  Applications.  North-* 
Holland  Publishing  Cg-." ,'  ~Afaf ft* V).  -  - 


form 


(3.11) 

£(hp)  ■  ha^  «AB 

where 

hd*f  (-det(hAB))1/2 

Substituting  (3.  20)  back  into  (3.  7)  then  leads  to 

(3.12)  p  WA  ♦  WA  a30  «gC  +  o^.g  ■  0. 

We  have  thus  established  the  following  result.  The  most  general  material 
momentum-energy  tensor  of  class  one  is  given  by 

(3.13) 


AB 

where  <7  is  a  symmetric  tensor  such  that 

(3.14)  v(oAB)  -a*®, 

(3.15)  bAB  VQ  ■  i*  VA 

A  AB 

if  and  only  if  •  0,  and  the  quantities  p,  W  ,  and  a  are  such  that 

AB 

(3.16)  £(h|i)  -  ho*°  «AB, 

(3. 17)  p  WA  +  WA  oBC  ^  +  oABb  -  0. 

4.  PHYSICAL  INTERPRETATIONS 


The  question  naturally  arises  as  to  Just  what  physical  interpretations 

A  AB 

can  be  attached  to  the  quantities  p,  W  ,  and  o  .We  take  here  the  view¬ 
point  that  physical  interpretations  should  be  consequences  of  the  mathematical 
model  employed. 

A 

We  first  show  that  the  vector  field  W  may  be  interpreted  as  a  velocity 

A  jo 

field.  Since  W  is  assumed  to  be  known  throughout  the  region  *  we  may 

define  a  system  of  fibers  in  ^  by 

A  AD  A 

W  is  uniquely  determined  by  T  and  W  W  ■  1. 


(4.1) 


dxA/dq  »  WA(x) 

where  t)  la  the  fiber  parameter.  Computing  the  differential  element  of  arc 
length  along  any  fiber  in  we  are  leada  to 
da2  •  hAB  WA  WB  dr\2  ■  dq2. 

Hence  (4. 1)  ahowa  that  WA  aatiaflea  the  definition  of  a  velocity  field  in  gen¬ 
eral  relativity. 

A 

The  next  queation  is  what  phyaical  quantity  admlta  W  aa  its  velocity 
field.  From  (2. 2)  we  obtain 

(4.2)  ti*TABWAWB. 

Transforming  to  a  frame  of  reference  in  /'  such  that* 

(4.3)  W  -5*, 

equation  (4. 2)  reduces  to  p  *  TQQ.  Tima,  since  (4.  2)  defines  a  reat  frame 

A 

relative  to  the  velocity  field  W  ,  and  the  T^q  appearing  above  is  the  energy 
in  this  rest  frame,  **  the  scalar  character  of  p  shows  that  p  is  the  rest-energy 
in  all  frames  of  reference.  Writing  (3. 16)  in  the  equivalent  form 

<*  hwA),A*h<’AB,AB 

A 

shows  that  p  h  W  is  a  flux  of  rest-energy  density  that  is  created  in  an  amount 

Aft  12  3  4 

equrfl  to  a  AB  per  unit  of  geometrical  volume  bdx  dx  dx  dx  .  From  the 

obvious  analogy  with  fluid  mechanics,  equations  (3.17)  may  be  Interpreted  as 

- * - 

A  This  can  always  be  done  for  any  vector  field  with  the  properties  exhibited 
by  W  :  see  E.  Goursat,  i .•mix*  «ur  u  paahlluaa  Herman  et  Cie., 

Paris  (1922),  p.  117.  The  symbol  ■  will  be  used  to  denote  evaluation  in  this 
framed  reference.  A  # 

From  (4. 3)  and  W.  W  »1  we  obtain  h  •  1.  Hence,  Tnn  in  the 
frame  defined  by  (4. 3)  is  numerically  equal  to  thr'NUakowBkian  value H>f 

T00~the  ««»rgy. 


AB 

describing  the  motion  of  a  fluid  with  mass  density  p  and  stress  tensor  o  , 

where  the  mass  density  is  created  in  conformity  with  equation  (3. 16).  The 
AB 

interpretation  of  o  as  a  generalized  stress  is  also  operationally  admissible. 

AB 

This  follows  from  the  fact  that,  by  (3. 14),  the  support  of  o  is  the  generic 

three-dimensional  space  U-  orthogonal  to  WA,  and  hence  oAB  may  be 

A  *  A 

determined  in  the  frame  W  *  6^  by  only  space-like  measurements.  Thus, 

since  the  quantity  p  has  been  shown  to  be  the  rest-energy,  we  may  interpret 

the  vector  field  WA  as  the  velocity  field  associated  with  the  flow  of  rest— 

y  AB 

energy  in  the  region  and  cr  as  the  generalized  stress  that  gives  rise  to 
this  energy  flow. 

The  above  interpretations  have  been  drawn  in  a  manner  which  is  in¬ 
dependent  of  what  actual  physical  processes  occur  in  Zd>  the  only  proviso  being 
AB 

that  T  be  a  material  momentum-energy  tensor  of  class  one.  We  may  thus 

* 1 

view  the  dynamical  processes  in  any  region  ^  with  a  material  momentum- 
energy  tensor  as  the  flow  of  rest-energy  regardless  of  what  actual  physical 
processes  occur.  The  continuum  in  relativity  theory  thus  turns  out  to  be  an 
energy  continuum— thus  the  name  relativistic  energy  mechanics. 

5.  BOUNDARY  CONDITIONS 

The  next  problem  to  be  resolved  in  the  resultant  energy  mechanics  is  that 

of  formulating  consistent  boundary  conditions .  Let  tZ  be  a  region  of  which 

supports  an  energy  mechanics  (i.e. ,  the  associated  momentum-energy  tensor 

of  Z?  is  material  and  of  class  one),  and  assume  that  exterior  to  TAB 

*The  generic  space  ~ZC  becomes  a  proper  space  if  and  only  if  »(W,  , )  ■  0: 

see  D.  G.  B.  Edelen,  "Rotation  Tensors  and  Irrotationai  Motions  in  '  J 
Einstein— Riemann  Spaces,"  Proc.  Nat.  Acad.  Scl.  (1963)  (to  be  published). 


is  not  material.  If  we  denote  the  boundary  of  2$  by  d<^  and  the  outward  normal 
to  9?3  by  NB,  then  a  necessary  condition*  for  the  solution  of  the  Einstein  field 
equations  is 

(5.1)  (TAB]N  «0,  <[TAB]dSf  -TAB). 

4*  —  B 

AD  AD  AD 

Here  T  (T  )  denotes  the  value  of  T  just  exterior  (interior)  to  82^  . 

"ab 

Under  the  above  assumptions,  we  may  use  (3. 13)  to  evaluate  T  .  When  this  is 
substituted  into  (6. 1),  we  obtain  the  conditions 
(5.  2)  ft  WA  W®  +  ff AB)  Nb  ■  T*®  Nfi. 

Hence,  if  we  operate  on  both  sides  of  (6. 2)  with  »,  set 
(5.  3)  FA  ■  r(TAB  N0), 

AB 

and  use  the  projective  invariance  of  a  as  stated  by  (3.14),  we  are  led  to 

(5.4)  o  ^  f(Ng)  •  FA. 

Equations  (6.4)  arc  just  the  relativistic  analog  of  the  classical  stress  boundary 

A 

conditions,  as  is  easily  seen  by  referring  (3.4)  to  a  rest  frame  relative  to  W 

\ 

Combining  (5.  2)  and  (5.4),  the  corresponding  boundary  conditions  of  the  W  are 
easily  obtained;  namely 

(5.5)  p  (W®  Ng)  WA  -  TAB  Ng  -  FA. 

Hence,  since  i*  A  0,  a  necessary  and  sufficient  condition  for  W  to  be  tangent 
to  8 S  is 

Na(1-*)<T  Ng)  -  0. 

The  evaluation  of  the  corresponding  conditions  for  the  h^g's  and  the  resulting 
* 

Edelen,  D.  G.  B.,  and  T.  Y.  Thomas,  Discontinuities  in  the  Einstein 
Field  for  General  Momentum-Energy  Tensors, "  Arch.  Rational  Mech.  Anal.  9 
(1962),  pp.  153-171. 


surface  dynamics  has  been  treated  elsewhere,  * 
6 .  CONSTITUTIVE  RELATIONS 


The  relativistic  energy  mechanics  established  in  the  previous  sections  lacks 

AB 

constitutive  relations  whereby  the  generalised  stresses  o  are  related  to  the 
velocity,  the  velocity  gradient,  and  any  other  fields  which  may  be  present.  Such 
relations  are  a  necessary  part  of  the  theory,  for  without  them  the  Einstein  field 
equations  and  the  basic  equations  of  relativistic  energy  mechanics  would  be  in¬ 
sufficient  to  determine  a  tenaorlal  solution  manifold.  * 

DEFINITION  3.  A  relation  of  the  form 

(6. 1)  o^  •  Q^O Vc,  WE.F,  h^,  . . .) 

is  said  to  be  a  constitutive  relation  in  the  relativistic  energy  mechanics. 

The  point  of  departure  in  constructing  such  relations  is  the  data  supplied 
AB 

by  section  3;  namely,  o  is  symmetric  and  such  that 

(6.2)  »(oAB)  ■  aAB, 

(6.3)  oAB  Vfi»  (iVA 

AB 

if  and  only  if  ■  0.  Since  (6. 2)  requires  the  support  of  o  to  be  the  generic 

three-dimensional  space  ^  orthogonal  to  the  vector  field  WA,  it  would  appear 

natural  to  confine  our  attention  to  quantities  whose  support  is  also  2^  . 

Edelen,  D.  G.  B.,  and  T.  Y.  Thomas,  "The  Dynamics  of  Discontinuity 
Surfaces  in  General  Relativity,"  (to  be  published). 

•• 

The  Einstein  field  equations  together  with  basic  equations  of  energy 
mechanics  constitute  a  system  of  11  independent  equations  for  the  determination 

A  AB 

of  the  variables  h._,  p,  W  and  o  among  which  there  are  only  23  independent 
a®  AB 

ones.  Thus,  if  we  append  ten  equations  for  the  determination  of  the  ten  a  *s 

without  introducing  any  new  variables  the  solution  manifold  would  be  tenaorlal 
because  of  the  occurrence  of  four  arbitrary  functions  of  position,  the  latter  being 
determined  by  the  choice  of  the  coordinate  system. 


Physically,  this  would  mean  that  an  observer  in  a  frame  of  reference  at  rest 

A 

with  respect  to  the  velocity  field  W  could  determine  both  constitutive  relations 
and  stresses  by  measurements  of  only  space— like  quantities.  This,  however,  is 
exactly  the  state  of  affairs  in  the  classical  mechanics  of  deformable  contlnuua. 

We  are  therefore  led  to  the  following  definition. 

DEFINITION  4.  A  relation  of  the  form 

(6.0  «AB-cAB<;t»  ,iB>.  ;<hAB).  -.) 

is  said  to  be  a  classically  admissible  constitutive  relation  in  the  relativistic 
energy  mechanics,  and  the  resulting  energy  mechanics  is  said  to  be  classically 
admissible. 

It  is  to  be  expected  that  only  in  the  cases  of  classically  admissible  con¬ 
stitutive  relations  will  the  relativistic  energy  mechanics  reduce  to  classical  con¬ 
tinuum  mechanics  by  an  appropriate  limiting  process. 

The  simplest  kinds  of  classical  continva  are  those  referred  to  as  isotropic. 
Although  the  idea  of  Isotropy  is  well  defined  in  the  classical  theory,  there  are 
several  possible  formulations  available  in  the  space-time  of  general  relativity. 

The  simplest  and  most  obvious  intuitively  is  that  associated  with  the  idea  of 
coincidence  of  the  principal  directions  of  stress  and  the  velocity  gradients.  A 
necessary  and  sufficient  condition  for  coincidence  of  the  principal  axes  of  o^ 

and  WA  0  is  that  crAB  be  a  power  series  in  W B  and  WQ  A,  with  scalar 

AB 

functions  of  W’A  B  as  coefficients.  This  would  violate  the  symmetry  of  a 
Hence  we  can  at  most  require  that  the  principal  axes  of  oAg  and  coincide. 

Applying  the  same  reasoning  in  conjunction  with  the  condition  (5. 2)  shows  that  we 


-12- 

can  at  moat  raquire  colncldanca  of  the  principal  axes  of  o._  and  »(W. .  ,). 

An  w  (A;B) 

Hanca  we  ara  lad  to  the  following  dafinltion  of  isotropy. 

DEFINITION  5.  A  ralatlviatic  energy  continuum  la  isotropic  if  and  only  if  tha 


axaa  of  a  Ag  coincide  with  tha  axes  of  «AB> 


If  wa  aat 


(8.5) 

where 


B  CB  AB  ,  CA  .  OB 

•a-'aC15  ••  *‘cDk  k  • 


k*8  .  «<hAB),  k._  • 


then  the  most  general  constitutive  relations  for  an  isotropic  continuum  which 
satisfy  the  conditions  ■  0  and  (8. 2)  are  given  by 


(87) 


where  the  coeffieiente  #A  are  scalar  functions  of  Since  detU^)  ■  0 

and  «AB  satisfies  its  own  characteristic  equation,  we  have 


(6.8) 


where 


(6.9) 


AC  D  B 


AC  B 


*D  +  !1  *  *c  +  *2  *  ~  l3  °* 


»  A 

h '  V 

.1-2  A  BCD  FGHK  .  . 

*2  4!  h  *  *  hAF  “BG  *CH  bK* 

.1-2  A  BCD  FGHK  . 

3  l!h  *  *  hAF  *BG  *CH  ^DK* 


A  BCD 


and  e  is  the  Levi-Civita  indicator.  Hence,  since  the  W'A  B  dependence 

of  the  scalar  functions  *,  can  be  written  in  terms  of  the  invariants 


(6.10) 


I  ^  *  *  ;A» 

I.,.1  ABCD  FGHK 
J  U2  H  h  hAFhBG'VC;H^D;K' 

i  »2  .  ±  h“2  eABCD  eFGHK  h  WWW 
j  3  4'.  *  *  “AF  WB;G  WC;H  WD;K* 

l“«  *  h"2  “•«*A1B»* 


/‘at  9*4. 
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the  moat  general  Isotropic  energy  continuum  admits  the  constitutive  relations 
(6.11)  a  •  jg  k  ♦  *  *C* 

Where 

<8  l2)  *A*fACl*  lv  V  ni»  n2»  U3‘  n4*  ‘  *  )*  A"°>1,2, 

and  the  f's  are  scalar  function*  of  their  indie  ated  arguments.  If  we  require  the 

AB 

continuum  to  be  claeelcally  admlaaible,  than  a  is  again  given  by  (6. 11),  but 
the  f  are  no  longer  functions  of  the  I's,  and  the  h'e  appearing  in  (6. 9)  are 
replaced  by  k's. 

We  still  have  to  satisfy  the  requirement  (6. 3).  Let  e  (i  •  1,  2,  3)  be 
defined  by 

AQ  A 

(6. 13)  *  a  V„  •  e  V  ,  (i  not  summed), 

iB  i  i 

where  v(VA)  *  V A.  Since  the  projective  invariance  of  the  V's  implies 

V.kAB-VB,  (6. 11)  gives 
iA  i 

(6.14)  oAB  V  ■  (m-  ♦  ®.  e  +  ?  e2)  VA  (i  not  summed). 

iB  *0  *1  i  2  i  i 

Hence  we  must  require  that  the  © 's  and  the  e  be  such  that 

i 

(6.15)  *0  +  1*1  *  ^2  *2  *  *  ’ 

It  the  above  results  are  examined  in  a  frame  of  reference  at  rest  with 
respect  to  the  velocity  field  WA,  the  classically  admissible  constitutive  re¬ 
lations  are  easily  seen  to  be  exactly  the  results  known  in  the  classical  theory. 

As  an  example,  consider  the  relation 

(6.16) 

where  k  and  q  are  constants.  This  is  clearly  classically  admissible.  An 
obvious  calculation  based  on  (6. 16)  yields 


AB  _  .  A  v  C.  .AB  AB 

O  ■  (p  +  k  «c)  k  +  2  q  <  , 
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AR  AR  A 

(1. 17)  °r  -a  -a  kAB  ■  3p +(2ij  ♦  SO  «*  . 

We  tins  obtain  an  immediate  baa  la  for  examining  the  energetic  analog  of  a  general 
viscous  fluid  if  we  make  the  etokes  aaaumptien  that  2  r\  ♦  3  V  >0. 

It  muat  be  clearly  borne  in  mind  that  we  are  dealing  with  an  energy  con¬ 
tinuum.  The  function  p  is  thus  not  necessarily  the  relativistic  analog  of 
hydrodynamic al  pressure.  If  p  were  assumed  to  be  the  pressure,  as  in  the 
usual  discussions,  certain  problems  arise  for  which  there  is  no  simple  answer. 
Einstein  clearly  saw  these  problems  and  chose  to  follow  an  interpretation  similar 
to  that  taken  in  this  paper.  "This  [referring  to  p)  must  not,  however,  be  con¬ 
fused  with  a  hydrodynamic  al  pressure,  as  it  serves  only  for  the  energetic  pre¬ 
sentations  of  the  dynamical  relations  inside  matter.  "* 

7.  STRESS  POTENTIALS 

Many  problems  in  classical  continuum  mechanics  are  significantly  simplified 
when  the  material  medium  is  such  that  it  admits  a  stress  potential  £.  The  con¬ 
dition  for  this  to  be  the  case  is 

for  all  dj,  where  tj  is  the  Cauchy  stress  tensor,  dj  is  the  rate  of  deformation 

tensor  (following  the  motion),  p  is  the  mass  density,  pQ  is  the  value  of  p  in 

the  reference  state,  and  £  is  the  time  derivative  of  £  (following  the  motion).  ** 

In  the  hope  of  achieving  a  similar  simplification  for  problems  in  relativistic  energy 

mechanics,  we  shall  look  for  a  suitable  generalisation  of  (7, 1). 

4 

Einstein,  A. ,  The  Meaning  of  Relativity.  Third  Edition,  Princeton  University 
Press  950),  p.  103. 

See  C.  Truesdell,  "The  Mechanical  Foundations  of  Elasticity  and  Fluid 
Mechanics. "  J.  Rational  Mech.  Anal..  1(1952),  p.  170. 


The  first  problem  to  be  resolved  is  thst  of  obtaining  an  acceptable  genera¬ 


lization  of  the  left-hand  side  of  (7. 1).  Since  (7. 1)  is  an  equation  following  the 
motion  of  the  classical  continuum  (Eulerian  specification),  we  refer  the  energy 

A 

continuum  to  a  frame  of  reference  at  rest  with  respect  to  the  velocity  field  W 

(that  is,  WA  ■  6q  ).  In  this  frame  of  reference,  x°  ie  the  natural  time  variable 

of  the  material.  Hence,  differentiation  with  respect  to  x®  is  the  relativistic 

analog  of  the  classical  time  differentiation  following  the  motion.  We  now  need  to 
0 

express  x  differentiation  in  the  rest  frame  by  means  of  an  invariant  differentiation 

A  dxA  0  * 

process.  Now  W  «  —  ;  hence,  in  the  rest  frame,  dx  •  ds.  We  thus  have 

. 

established  the  correspondence  E  -  -jjj-  ,  for  A  a  scalar  function.  In 

classical  mechanics,  the  density  p  is  the  reciprocal  of  the  element  of  converted 

volume.  Thus,  if  denotes  the  slice  of  $  over  which  the  initial  data  is  specified, 

and  hQ  denotes  the  image  of  h  in  >'■  under  the  inverse  of  the  motion  defined  by 

A 

the  trajectories  of  the  W  field,  we  have  the  natural  correspondence 
p  p  1  ■*  hQ  h  Combining  the  above  results,  we  obtain 
(7.J)  ,,-'t 

Because  of  the  occurrence  of  the  ratio  hQ  h  *,  it  would  seem  appropriate  to  refer 
to  A  as  the  convective  stress  potential  (classically,  the  quantity  E  should  also 
be  referred  to  as  the  convective  stress  potsntial  due  to  the  factor  p  p~V 
Turning  to  the  right-hand  side  of  (7. 1),  we  note  that* 

(7.3)  2  *<A;B)  * 

and  hence 
- • - 

See  J.  A.  Schouten,  Riccl-C aiculus  (Second  Edition),  Springer— Verlag, 

Berlin  (1954),  pp.  103-111. 
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(7.4>  2  «AB  *  w  (h^). 

Evaluating  (7.4)  In  the  rest  frame  given* 

(7' 5*  2  *AB  “  <hAB  “  hA0  bB0),  0‘ 

#  * 

(Remember  that  hQ0  •  1 . )  The  teneor  «A0  le  thus  a  time  derivative  following 
the  motion.  The  question  is  whether  h^  —  hQA  hQ0  may  be  looked  upon  ae  an 

A 

Eulerianrate-of-etrain tensor  in  the  rest  frame.  Since  W  is  a  time-Uke  unit 

vector,  we  may  pick  a  system  of  four  vectors  v  (b  ■  1,  . . . ,  4)  and  a  reciprical 

b  b 

system  vA  such  that 

A  *  «A  A  *  .a  .  A  „  A 
(7.8)  v  vb-6b,  v  vA-6b,  and  v  . 


We  then  have 

(7. 7)  h 


0  0  3  i  i 

AB  ‘  “o  VA  ’B  +  *  \  VA  VB* 

A..B 


Using  (7.6)  and  the  fact  that  hAB  W  W  *1  gives  0Q  •  1,  and  hence  we  have 

*  3  i  i 

(hAB-hA0  W‘  *1  ’A  V 


i-i 


Thus,  noting  that  w  (vA>  ■  vA  and  that  the  are  the  nonholomonic  components 


A  i  i 

of  h^  in  the  generic  space  normal  to  W  ,  we  see  that  vA  v0  may  be 

looked  upon  as  a  classical  strain  tensor  in  the  three-dimensional  space  the 

exception  being  that  the  vA  are  in  general  nonlntegrable  (l.e. ,  vA  ¥  u 1  A>.  The 

tensor  2  «A0  is  thus  seen  to  be  the  natural  generalisation  of  dj. 

In  view  of  the  above  results,  we  have  the  following  characterization  of  a 


convective  stress  potential  in  general  relativity. 

*Edelen,  D.  G.  B.  and  T.  Y.  Thomas,  "On  the  Characterisation  of  Born 
and  Contemporaneous  Rigidity  and  the  Question  of  Their  Equivalence, "  Tensor 
(to  be  published). 
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DEFINITION  6.  A  scalar  function  A  la  said  to  be  a  convective  stress  potential 
of  a  relativistic  continuum  if  and  only'  if 

““  Jf  •  2  “AB  -AB 

!■  an  Identity  in 

The  explicit  occurroact  of  the  factor  hQ  in  (7. 8)  render*  this  equation 
difficult  to  deal  with.  For  this  reason,  w*  transform  (7.8)  into  an  equivalent 
form  which  does  not  contain  hQ  explicitly.  One  may  easily  verify  that 

(7.9)  £(hexp(—  «  ds))  ■  0, 
where 

(7.10)  <  * -AB  ^  ‘ 

A 

and  the  integration  is  to  be  performed  along  the  trajectories  of  the  W  field. 

An  integration  of  (7. 9)  thus  leads  to 

(7. 11)  h  exp(—  j  «  ds)  •  hQ. 

In  view  of  (7. 11),  we  introduce  the  scalar  density  P  by 
def 

(7.12)  P  •  h  A  exp(—  <  ds). 

v  / 

Now  •£  (A),  and  hence,  by  (7.9)  and  (7. 11), 

£  (P)  •  h  exp<-  j  «  ds)  £(A)  ■  hQ  ~  . 

Substituting  from  (7.8)  into  the  right-band  side  of  the  above  equation,  we  establish 
the  following  result. 

A  convective  stress  potential  A  exists  if  and  only  if  there  exists  a  scalar 
density  P  such  that 

(7.13)  £  (P)  •  2h  <yAB 


-18- 


1«  an  identity  in  If  this  condition  is  satisfied,  than  A  is  given  by 

(7. 14)  A  ■  h”1  P  np( J  «  ds). 

For  obvious  reasons,  ws  shall  refer  to  P  as  the  stress  potential. 

8.  STRESS  POTENTIALS  AND  CONSTITUTIVE  RELATIONS— AN  EXAMPLE 


It  Is  shown  In  the  Appendix*  that  necessary  and  sufficient  conditions  for 
the  existence  of  a  stress  potential  is  that  the  following  equations  be  intsgrsble 


for  P: 

(8.1) 

h  oAB  ■  »(P,.  ), 

hAB 

(8.2) 

K  +  2  AA  WA  ■  0. 

A 

AH 

If  P  exists  (that  is,  if  a  satisfies  the  conditions  (A.  IS)),  equations  (8. 1)  are 
a  system  of  10  constitutive  relations  which  determine  the  generalized  stresses. 
Now,  if  P  depends  only  on  bAB  and  «AB,  the  generalized  stresses  given  by 
(8. 1)  will  be  uniquely  determined  functions  of  the  dynamical  state  of  /S '•  In  this 
case  (8.  2)  represents  an  additional  condition  to  the  Einstein  theory.  On  the  other 
hand,  If  P  depends  on  at  least  one  other  function,  the  relation  (8.  2)  Is  consistent 
with  the  Einstein  theory;  in  fact,  it  provides  a  natural  additional  constitutive 
relation. 

We  shall  illustrate  this  result  by  the  following  example.  Let 

(8.3)  P  ■  —  2h  (p  —  X  «); 
then  (8. 1)  and  (8. 2)  give 

(8.4) 
and 


AB  .AB,  v.  _ ,  CA  BD 

o  ■  -  k  (p  -  X«)  -  2  X  «CD  k  k 


0  -  -  £  (p)  ♦  «  £  (X)  +  XhCD  £  («CD)  , 


Apply  the  substitutions  E  ■  P  and  •  ■  1,  to  (A.  2)  to  obtain  (7. 13). 


(8.5) 
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respectivoly  now,  we  have 

£(«)  -  hCD  £  (4CD)  4  «CD  £  (hCD), 
so  that  (t.  5)  may  b«  writtsn  In  ths  equivalent  form 

(8.6)  £(p-X«)  -XhAE  hEF  *AB  t  (hgj,). 

Also,  however,  we  have 

2*EF  "  f  (L<hEF0» 

and  hence  a  direct  calculation  leads  to 

.AE.BF  „  AE.BF 

h  h  *AB  £  ^EF  2h  h  ‘AB*EF* 

Thus  (8.5)  Is  equivalent  to 

(8.7)  £  (p-X«)  -  2XhAEhBF*AB«Ep. 

Since  p  —  U  is  a  scalar,  Integration  of  (8.  7)  leads  to 

<8. 8)  P  "  P0  “  Vo  +  +  2  j  KhAEhBF'AB‘EFd>» 

A 

where  the  integration  is  to  be  performed  along  the  trajectory  of  the  IV  field 
from  the  appropriate  point  in to  the  point  in  &  at  which  the  left-hand  side 
of  (8. 8)  is  to  be  evaluated.  (Here,  pQ,  XQ,  and  cQ  are  the  values  of  p,  \,  and 
<  in  JnSt  Hence,  for  given  X,  both  the  generalised  stresses  oAB  and  the 
function  p  are  uniquely  determined  functions  of  the  dynamical  state  of  . 

On  the  other  hand,  if  we  were  to  set  p  •  0  and  consider  X  as  given,  the  above 


(8. 10) 
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In  this  case  we  again  obtain  a  determination  of  the  generalized  stresses,  but  only 
for  those  motions  which  are  such  that  (8. 10)  is  satisfied. 

Consider  now  the  case  for  which  X  •  0.  This  gives 

(8.11)  P  •  -  2  h  p,  crAB«-pkAB, 
and 

(8.12)  £  (p)  ■  0. 

AB 

With  cr  given  by  (8. 11),  the  basic  equation  (3. 16)  leads  to 

(8.13)  £(h  ji)  ■  —  p  £  (h). 

Thus,  applying  the  result  (8. 12)  we  have  the  following  conclusion.  The  existence 

AB 

of  a  stress  potential  for  -p  k  implies  the  conservation  law 

(8.14)  (p  WA).a  *  0,  p  ■  n  +  p. 

Combining  this  conclusion  with  the  results  of  Thomas*  we  see  that  a  region 
of  £  with  the  material  momentum-energy  tensor 

A  R  AB 

tab-'w  w  “Ph  *  +p>  0 

AB 

contains  a  material  point  whose  trajectory  is  a  geodesic  if  the  stress  tensor  -p  k 
admits  a  stress  potential. 

T.  Y.  Thomas,  On  the  Geodesic  Hypothesis  in  the  Theory  of  Gravitation,  " 
Proc.  Nat.  Acad.  Sci.  (U.S.A.),  48  (1982),  pp.  1567-1569. 
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9.  PATH  DENSITY 

The  statement  of  a  material  momentum-energy  tensor  of  class  one  is 

AB 

seen  to  require  the  statement  of  the  generalized  stress  tensor  c  .  In  general, 

the  stress  tensor  characterizes  the  resistance  of  the  energy  continuum  $  to 

deformations;  as  such,  this  tensor  depends  on  the  particular  properties  of  the 

energy  continuum  under  consideration.  Statement  of  cAB  is  thus  not  an  appended 

condition  to  the  Einstein  theory;  it  serves  to  characterize  the  actual  properties  of 

the  particular  energy  continuum  under  investigation. 

We  shall  henceforth  assume  that  the  cAB  are  known  functions  of  the 

state  of  me  on  of  $  .  Since  d^B  has  the  generic  three-dimensional  space  4^ 

as  support,  this  tensor  will,  in  general,  depend  on  k^,  »A,  k  ,  and 

as  well  as  any  other  fields  which  are  considered  in  the  particular  problem  under 

AB 

investigation.  Hence,  under  solution  of  the  Einstein  equations,  the  tensors  c 

and  <  .  _  will  be  known  functions  along  the  trajectories  of  the  wA  field. 

AB 

The  state  of  motion  of  an  energy  continuum  &  is  determined  by  solving 
the  Einstein  field  equations  together  with  the  basic  equations  of  energy  mechanics 
given  in  Sec.  3.  In  order  that  this  may  be  accomplished,  appropriate  initial 
(Cauchy)  data  must  be  specified.  Let  7/ be  the  section  of  $  over  which  the  initial 
data  is  specified.  In  addition,  let  5  be  the  assigned  value  of  u  on  */n  . 

A  function  fundamental  in  our  consideration  is  the  path  density  S.  This 
function  is  a  scalar  density  defined  on  the  trajectories  of  the  field  by  the 

conditions 

•The  generic  space  V  is  unique  owing  to  the  fact  that  the  vector  field 
is  unique. 
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(9.0  ns)  -  h  oAB  «AB,  h  ■  (-dct(hAB»1/2, 

and 

(9.2)  Sj  /  0. 

l/n# 

By  definition  of  the  Lie  derivative  of  a  scalar  density,  we  have 

«9-3'  £(S>.4StswA  fUs.. 

ds  ;A  ds  ’ 

where  s  is  the  arc-length  parameter  along  the  trajectories  of  the  field 
and 

(9.4)  « 

Substituting  (9.  3)  into  (9.  1)  and  using  the  condition  (9.  2)  to  evaluate  the  arbitrary 
function  which  arises  in  the  integration  process,  we  obtain  the  following  explicit 
evaluation  of  the  path  density: 

(9.5)  S  ■  exp(-  f€ds)  j exp(  ds)  h  cAB  « Afi  ds, 

the  integrations  being  performed  along  the  trajectories  of  the  WA  field  from  the 
appropriate  point  in  to  the  point  under  consideration.  Since  the  generalised 

stresses  are  assumed  to  be  known  on  the  trajectories  of  the  field,  (9.  5) 
determines  a  unique  function  on  the  trajectories  of  . 

There  are  several  transformations  of  the  path  density  S  which  will  be 
useful.  If  we  set 

(9.6)  S  *  hQ, 

then  Q  is  a  scalar  function  and  (9. 1)  becomes 

(9.7)  £  (Q)  •  « AB . 

Proceeding  in  the  same  manner  as  above,  we  then  obtain 
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(9.8)  Q  ■  expt-Jtds)  fexp(J«d8)  cAB  «AB  ds  . 

The  function  Q  will  be  referred  to  as  the  path  function.  If  we  set 

(9.9)  S  ■  qhexpf-  J^ds) , 

then  q  is  a  scalar  function.  By  the  Identity 

(9.10)  £(hexp(-J«ds))  b  0, 
equation  (9.  1)  becomes 

(9.11)  exp(-J« ds) £ (q)  * 

Integrating  (9.  11),  we  then  obtain 

(9.12)  q  ■  Jexp(  r«ds)  aAB  «AB  ds. 

This  function  will  be  referred  to  as  the  path  kernel. 

1 0 .  CONSERVED  Q  UANTITIES 

Knowledge  of  the  path  density  S  allows  us  to  prove  the  following 
fundamental  result. 

The  Einstein  field  equations  imply  the  existence  and  uniqueness  of  the 
scalar  density 

,  »  def 

(10. 1)  t  »  hu-S 

such  that 

(10.2)  £(♦)  !  (IW*)  *  0. 

» A 

The  proof  is  as  follows.  Let  S  be  given  by  (9.  5);  then  (9.  1)  holds. 
By  (3.  18),  the  Einstein  field  equations  imply 

(10.3)  £(h  u)  -  «AB«ab- 
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AB 

Eliminating  the  expression  e  between  (9.  0  and  (10.  3)  gives 

(10.4)  £  ( u  h  -  S)  *  0. 

Now,  u  is  an  eigen  value  of  T*B,  so  that  k  u  +  -■  R  is  an  eigen  value  of  the 
Ricci  Tenson.  *  Hence  a  is  a  differential  Invariant  of  the  Einstein— Rlemann 
space  under  consideration.  The  density  uh  is  thus  uniquely  determined  by  the 
geometry  of  the  space-time  continuum.  The  path  density  S  is  also  unique. 

Hence  the  scalar  density  #  ■  uh  -  S  is  unique  and  satisfies  (10.  2). 

There  are  several  transformations  of  4  which  will  be  useful.  If  we  set 

♦  *  ha,  then  x  is  a  scalar  function  such  that 

(10.  5)  0  *  (cpW^)  A  ■  ~  +  «  cp  ■  £(cp)  +  t  <p 

;a  <38 

and  cp  ■  u-Q.  In  obtaining  the  last  result  we  have,  of  course,  used  (9.  6).  If 
we  set  4  »  thexpl-j*  «ds),  then  4  is  a  scalar  function.  By  use  of  the  identity 

J  i 

(9.  10),  equation  (10.  2)  reduces  to  0  *  £(♦)  *  —  .  Hence  ♦  is  always  a  constant 
of  the  motion.  On  the  other  hand,  a  is  a  constant  of  the  motion  only  if  the  motion 
is  incompressible,  that  is  <  *  0 .  In  an  obvious  manner,  we  also  obtain 

*  ■  uexp(-J*«ds)  -  q,  when  use  is  made  of  (9. 9). 

1 1.  THE  EXISTENCE  AND  UNIQUENESS  OF  INTRINSIC  ENERGY 

As  a  consequence  of  Einstein's  fundamental  discovery  that  mass  and  energy 
are  equivalent,  any  change  in  the  energy  of  a  body  results  in  a  change  in  its  mass. 
Since  the  energy  of  a  body  is  a  function  of  its  dynamical  state,  so  is  its  mass. 
Fundamental  to  the  idea  of  a  material  body,  however,  is  the  existence  of  an  intrinsic 
mass  (or  energy)  which  is  Immutable.  Accordingly,  when  discussing  material 

— *This  1*  an  immediate  consequence  of  the  Einstein  field  equation 
RAB~RhAB^  "  ^AB' 
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bodies  in  the  general  theory  of  relativity,  additional  conditions  are  often  imposed 
so  as  to  secure  an  intrinsic  mass.  This  section  shows  that  it  is  unnecessary  to 
impose  such  additional  conditions  -  the  Einstein  field  equations  imply  the  exist¬ 
ence  and  uniqueness  of  just  the  required  intrinsic  quantity. 

We  have  shown  that  if  a  region  V  of  an  Einstein-Rlemann  space  is 
occupied  by  a  material  body  (that  is,  the  associated  momentum-energy  tensor  of 
G  is  material  and  of  class  one),  then  there  exists  a  unique  scalar  density  ♦  on 
G  such  that 

01.  1)  <*  W^)  A  »  0 

and 

(11.2)  ha  *  i  +  S. 

Equation  (11.  1)  states  that  the  flux  of  i  is  conserved  during  the  motion  of  the 
body.  Since  h  u  is  the  density  of  rest-energy  (that  is,  the  energy  which  would 
be  observed  in  a  coordinate  system  such  that  WA  *  6^),  equation  (11.2)  states 
that  the  density  of  rest  energy  is  uniquely  representable  as  the  sum  of  a  path- 
dependent  density  S  (the  path  density)  and  a  density  i  which  is  conserved.  Now, 
by  (9.2)  and  (11.2),  we  have 

Thus,  *  is  a  density  of  energy.  Since  t  la  conserved,  it  is  intrinsically 
associated  with  the  body.  Hence,  the  Einstein  theory  implies  the  existence  and 
uniqueness  of  an  intrinsic  energy  density  of  a  material  body,  namely  i. 

The  functions  cs  and  *  introduced  in  the  last  section  may  now  be 


*  Interpreted.  Since 
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<U.„  .-u-Q. 

and 

(11.4)  <ttwA)  A  *  °* 

we  may  Identify  cp  with  the  Intrinsic  reat-energy .  Noting  that  u  is  the  total 
rest-energy,  we  see  that  Q  is  the  expendable  rest-energy;  namely,  that  pari 
of  the  rest-energy  which  is  a  consequence  of  the  dynamic  processes  to  which 
the  body  is  subject.  From  the  identify 

(11.5)  £(h  exp(-  f«ds))  i  0 , 

we  see  that  exp(  J *  ds)  is  proportional  to  the  convected  element  of  geometrical 
volume;  that  is,  we  have  exp(  J(ds)  •  h/h^  where  h^  is  the  value  of  h  at 
the  corresponding  point  in  ^ .  Thus  since  $  »  ^exp(-J«  ds)  -  q  and  *  0, 
we  see  that  *  is  the  intrinsic  rest-energy  per  unit  of  convected  geometrical 
volume.  In  particular,  we  have  hg^uh-Shg)  *  q,  so  that  q  is  the  total 
change  in  rest-energy  density  per  unit  of  initial  geometrical  volume. 

12.  THE  QUESTION  OF  PATH  DEPENDENCE 

We  have  seen  that  the  path  density  S  is  given  by 
(12.1)  S  ■  exp(-J*«ds)  Jexp<  J«ds)  h  oAB  «AB  ds, 

and  that  this  quantity  together  with  h  u  serves  to  determine  the  intrinsic  energy 
of  an  energy  continuum.  In  general,  the  path  density  will  depend  directly  on  the 
trajectories  of  the  WA  field  in  addition  to  the  quantities  appearing  in  the  integrands 
of  (12.  1).  In  this  section  we  examine  the  conditions  under  which  S  is  path  inde¬ 
pendent  in  an  appropriate  sense. 

If  S  is  to  be  path  Independent,  we  must  require  the  integrand 
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Is  A3 

exp(j«da)ho  <  A0  to  be  a  total  differential.  We  thus  suppose  that  the  gener¬ 
alized  stress  field  is  such  that  there  exists  a  scalar  function  A  on  £  with  the 
property  that 

(12.  2)  ~  (h  exp(  J*«ds)  A)  ■  h  exp(  J*«  ds)  cAB  tA 

holds  for  all  kinematically  possible  motions  in  /&  .  If  this  supposition  is  satisfied, 
(12.  1)  gives 

(12.3)  S  *  h  A  —  exp(  —  J  «  ds)  hQ  AQ  , 

f 

where  h^  and  are  the  values  of  h  and  A  at  the  corresponding  point  in  *  h  &  . 

We  thus  have  to  examine  the  conditions  under  which  (12.  2)  holds  for  all 
kinematically  possible  motions  in  .  Now, 


^  [h  exp(  J €  ds)  A^  i  h  exp  ( J*  ds)  [«  A  +  ^  j  , 

and  hence  (12.  2)  is  equivalent  to 

(12.4)  oAB  dA  +  <A. 

AB  ds 

Noting  that  £(h)  *  «  h ,  we  finally  obtain  the  condition 

(12.5)  £(Ah)  *  hcAB  .  A  _ 

AB 

for  all  kinematically  admissible  « A0 .  This*  however,  is  exactly  the  requirement 
for  the  existence  of  a  stress  potential  P.  Hence,  setting  P  *  2hA,  we  have  the 
following  result:  A  necessary  and  sufficient  condition  for  S  to  be  path  Independent 
in  the  above  sense  is  that  there  exists  a  stress  potential  P.  If  this  condition  is 
satisfied,  then 

(12.6)  S  «  P/2 -exp(-;«ds)  Pq/2. 

We  have  seen  that  in  the  case  of  a  path  density,  we  have  the  intrinsic 


density  function 
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(12.7)  i  -  hu-S. 

Under  the  above  result,  (12.  6)  and  (12.  7)  give 

(12.8)  *  »  hu-P/2+  expt-J^ds)  PQ/2 

However,  Hexp(- j^ds)  PQ/2)*  0,  and  hence  the  quantity  hu-  P/2  is  an 
intrinsic  quantity.  In  addition,  since  P/2  is  path  independent  and  hence  a 
function  of  position  only,  we  have  the  local  conservation  law 

(12.9)  ((hu-  P/2)WA}  .  -  0 

f  A 

at  all  points  in  $ .  We  thus  see  that  path  independence  of  the  path  function 
implies  the  existence  of  a  local  conservation  law. 

13.  A  FUNDAMENTAL  DIFFERENTIAL  ON  THE  TRAJECTORIES  OF 
(THE  LAW  OF  THERMODYNAMICS?! 

As  a  consequence  of  the  Einstein  theory  and  the  requirement  that  a 

region  (3  have  an  associated  momentum-energy  tensor  which  is  material  and  of 

class  one,  the  dynamical  processes  interior  to  /3  can  always 

be  viewed  as  an  energy  mechanics.  In  the  customary  approach,  any  basic  energy 

mechanics  is  governed  by  the  laws  of  thermodynamics.  One  is  thus  led  to  inquire 

as  to  the  relations  between  the  previous  results  and  the  laws  of  thermodynamics. 

Although  it  may  have  seemed  disguised,  we  have  already  derived  the 

first  law.  One  of  the  basic  equations  of  relativistic  energy  mechanics  reads 

(13.  1)  £(hj  «  h  cAB  «AQ  . 

A 

Since  the  Lie  derivative  formed  from  W  is  the  relativistic  analog  of  the  time 
derivative  following  the  motion,  equation  (13.  1)  states  that  the  time  derivative  of 
the  rest-energy  density  following  the  motion  is  equal  to  the  power  developed  by  the 
generalized  stresses  per  unit  of  geometrical  volume  dx'dx^dx^x4  -  the  first  law. 
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In  order  to  proceed  further,  we  must  use  the  results  of  the  Appendix  to 
decompose  the  generalized  stress  tensor.  Set 

(13.2)  oAB  .  lAB  +  M^8  +  , 

AB  Ag  ig 

where  L  and  are  those  parts  of  o  for  which 

(13.  3)  f(P)  *  2  h  LAB  i  Ab 

v  -• 

(13.4)  £(£)*  2  6  h  N^8  « AB 

are  identities  in  «Ag»  and  I1^8  is  what  is  left.  Thus,  LAB  satisfies 
eABCD(l)  „  an(J  ^jAB  gati8fie8  gABCD^j  MFG  +  ECDFG(M)  ^B  + 

EFGAB^j  *  0  where  EA8^'E  (.  )  js  the  operator  defined  by  (A.  12). 

Now,  if  we  substitute  (13.  2)  into  the  right-hand  side  of  (13.  1)  and  use  (13.  4), 
we  are  led  to  the  expression 

(13.5)  £(hu)  «  £(P/2)  +  0-1  £(£/2)  +  * 
where 

(13.6)  *  »  h  l^8  « AB  . 

By  definition,  D(B)  s  £  (B)  ds  is  the  Lie  differential  of  any  geometric  object  B 

A 

with  respect  to  the  trajectories  of  the  W  field,  and  hence  is  a  differential 
following  the  motion  -  a  convective  differential.  Thus,  equation  ( 13.  5)  leads  to 

A 

the  following  fundamental  differential  on  the  trajectories  of  W 

(13.7)  D(hu)  *  D(P/2)  +  6_1  EKE/2)  +  tds  . 

Associated  with  this  we  have  the  intrinsic  quantity 

(13.  8)  1  *  hu-  P/2  +  exp(  -  ds)  Pq/2  -  ’  <  ds 

-|exp(-J«ds)  f  6  1  d(£exp(  <d*))  . 

We  now  have  the  desired  result.  Equation  (13.  7)  states  that  the  convective 
differential  of  the  density  of  rest-energy  is  equal  to  the  convective  differential  of 
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the  recoverable  energy  stored  In  the  deformations  of  the  energy  continuum  (i.  e. , 
the  energy  represented  by  the  path  Independent  potential  P/2)  plus  the  irrecoverable 
energy  increments  represented  by  the  convective  differential  form  6  1  D(Ll 2)  and 
yds.  The  convective  differential  form  B  1D(E/2)  is,  however,  exactly  analogous 
to  that  encountered  in  classical  thermodynamics  -  the  convective  differential 
operator  D  being  that  which  is  used  in  proper  formulations  of  thermodynamics. 

We  are  thus  led  in  a  natural  fashion  to  define  p  1  as  the  intrinsic  temperature  and 
E/2  as  the  intrinsic  entropy  density  of  the  energy  continuum.  In  fact,  if  we  insist 
that 

(13.  9)  D(E/2)  >  0  , 

then  (13.  7)  and  (13.  8)  state  the  relativistic  analogs  of  the  first  and  second  laws. 

A  general  material  region  of  an  Elnstein—Riemann  space  will  not 
necessarily  be  such  that  (13.  9)  is  satisfied.  The  relativistic  energy  mechanics 
is  thus  a  richer  subject  than  that  which  would  be  obtained  through  a  purely  thermo¬ 
dynamic  description.  This  added  richness  is,  however,  no  blessing  in  disguise  - 
if  anything,  relativistic  mechanics  is  ail  too  rich  for  our  analytical  blood.  We 
therefore  take  the  opportunity  afforded  by  the  above  result  to  constrain  the  universe 
of  discourse  of  our  mechanics. 

DEFINITION  7.  A  material  region  of  an  Einstein-Riemann  space  is  said  to  be 
thermodynamically  admissible  if  and  only  if  D(E/2)  >  0  . 

The  differential  form  n  ds  may  cause  some  pause  to  those  accustomed 
to  the  classical  notations  of  thermodynamics,  to  say  nothing  of  the  fact  that  D(hu) 
is  generally  not  an  exact  form.  The  inexactness  of  the  form  D(hu)  is  easily 
dispensed  with  on  noting  that,  contrary  to  classical  notions,  the  mass  of  a 
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relatlvistic  body  depends  intrinsically  on  Its  history,  and  hence  the  system  is 
not  locally  closed  with  respect  to  the  measure  u,.  If  one  desires,  the  existence 
of  intrinsic  energy  can  be  used  to  reformulate  the  description  in  terms  of  a 
system  which  is  closed  with  respect  to  its  history.  At  this  point  in  our  work 
we  can  only  conjecture  the  significance  of  the  form  »  ds.  If  1  *  0,  (13.  7) 

«.]  AO 

gives  D(hu)  =  D(P/2)  +  0  D(C/2).  Now,  D(hu)  *  ha  «ABd®  by  (3.  16), 
so  that  the  above  equation  becomes 

AT*  —1 

(13.10)  D(P/2)  *  he  «ABds-6  D(E/2)  . 

AB 

In  this  equation,  D(P/2)  is  an  exact  form  while  o  «Agds  has  the  dimensions 

,  ; 

of  work.  Hence  (13.  10)  gives  a  description  of  (o>  which  is  locally  closed  with 
respect  to  P/2  and  which  has  the  familiar  form  "dW  —  T  dS"  .  If,  however,  a 
system  emits  gravitational  waves,  it  should  be  Impossible  to  obtain  a  locally 
closed  system  with  just  the  usual  thermodynamic  variables  Hence  we  conjecture 
that  the  differential  form  t  ds  describes  the  thermodynamic  effects  of  gravitational 
radiation. 

14.  AN  EXAMPLE 

The  definitions  of  intrinsic  temperature  and  intrinsic  entropy  density 
given  in  the  last  section  are  of  a  formal  nature.  In  order  to  clothe  this  formality 
with  physical  significance,  we  examine  these  concepts  in  the  context  of  a  specific 
problem. 

Consider  a  material  region  $  of  an  Einstein— Riemann  space  with 
generalized  stress  tensor  given  by 
(14.  1)  a  «  -  pk  , 
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where  p  is  a  scalar  fvmctlon.  It  has  been  pointed  out  that  the  stress  tensor 
(14.  1)  admits  a  stress  potential  if  and  only  if  £(p)  s  ■  0  .  If  we  substitute 
(14.  1)  into  (3.  16),  we  have 

(14.2)  "  h-1  £(hp)d«f  x» 

where  p  ■  u  +  p  and  h  p  is  what  is  usually  identified  with  the  mass  density 
of  the  material  region. 

If  the  conclusions  of  the  last  sections  are  correct,  we  should  be  able 
to  obtain  all  significant  thermodynamical  variables  as  consequences  of  the 
Einstein  theory.  Let  'p  be  defined  by  'p  *  J  x  ds  ,  where  x  is  given  by 
h  ^lOip)  and  the  integration  is  to  be  performed  along  the  trajectories  of  the 
field.  We  then  define  the  scalar  "p  by  "p  *  p  -  'p  ,  so  that  p  »  'p  +  "p 
and  *  0.  The  decomposition  (13.  2)  gives 

(14.3)  oAB  «  -d’p  +  "p)  kAB  , 

from  which  it  may  easily  be  verified  that  the  equation 

(14.4)  £(P)  -  -2  hC'pkAB)  «AB 


is  satisfied  for  all  kinematically  admissible  with  P  *  -  2  h  "p  .  Hence, 

by  (13.  3)  and  the  above  result,  we  have  LAB  »  -  "p  kAB  and  MAB  *  —  'p  kAB  , 
•  q  ,  it  being  assumed  that  ♦  «  0  .  We  therefore  look  for  functions  6  and 


<L !£.  . 


ds 

E  such  that  (13.  4)  is  satisfied.  This  leads  to  the  relation 
(14.  5) 

Since  £(h)  ■  h«  ,  the  most  general  integral  of  (14.  5)  is 
(14.6)  E  »-2he'p+U, 


£(E)  -  2  h  6  m^ab  »  -  2  h  6  'p  « 


where 


(14.7) 


2  h  £('p  6)  ■  £(U)  ; 
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the  latter  equation  being  nothing  more  than  the  condition  stated  by  (A.  8). 
Equations  (14.  6)  and  (14.  7)  are  two  relations  between  the  three  unknowns  6,  E  , 
and  U.  We  must  therefore  have  one  additional  relation.  Now,  as  previously 
noted,  the  stress  tensor  (14. 1)  requires  the  specification  of  one  additional 
constitutive  relation  in  view  of  the  occurrence  of  the  function  p.  This  is 
exactly  reflected  by  the  need  for  one  additional  relation  for  the  determination 
of  E  ,  and  U.  We  thus  see  that  we  may  prescribe  the  function  U ,  and 
thereby  prescribe  the  required  constitutive  relation.  This  amounts  to  pre¬ 
scribing  the  equation  of  state  for  the  material.  In  the  simplest  case,  namely 
U  *  0,  (14.  7)  leads  to  p  n  «  pq  °o  ’  w^ch  may  be  written  in  the  more 


£  , 

T  v 


suggestive  form 

(14.8) 

where  the  definition  of  Intrinsic  temperature  is  used:  T  »  ^  *  intrinsic 
temperature.  If  we  take 

(14.9)  £<f>  *  Xh  d8>  £(P)  ~  'P  «  j  , 

where  ■  0  ,  then  (14.  7)  yields  pS  «  X  exp(-J«  ds)  .  Rewriting  this  with 
T  Instead  of  8  1  and  i .multiplying  by  h  gives  »  X  h  exp(-  J*  ds)  ,  so  that  we 
have  tty)  «  0  ,  the  state  equation  for  a  "perfect  fluid.  "  From  the  above 
examples,  we  see  that  not  only  do  the  ideas  introduced  in  the  last  section  make 
sense  for  the  case  of  a  general  gas;  they  also  provide  a  natural  manner  of  stating 
relativistically  valid  equations  of  state. 

As  a  last  note,  we  examine  the  Implication  of  the  requirement  that  the 
region  be  thermodynamically  admissible.  From  (14.  5)  we  have  £(E)  ■  —  2  h  8  'p  «  , 
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and  hence  the  requirement  (13.9)  reads  h  8  'p  «  <  0.  Now,  h  is  strictly 
positive  and  6  may  be  taken  as  strictly  positive  since  the  reference  value  of 
6  is  undetermined  by  (14.  7).  Hence,  (14. 1)  describe*  a  thermodynamically 
admissible  region  if  and  only  if  'p  «  <  0  ;  that  is,  if  and  only  if  a  positive 
"pressure"  can  only  produce  a  negative  dilatation. 


.  AB.  AB 
»(cr  )  *  C 


APPENDIX.  EXISTENCE  REQUIREMENTS 

AB 

Let  the  symmetric  tensor  a  be  given  and  such  that 

(A.  1) 

We  examine  under  what  conditions  a  scalar  6  and  a  scalar  density  £  exist 
such  that 
(A.  2) 


£(E)  «  2  6  h  oAB« 


*  AB  *  f(W(A;B))  r(£  hAB)/Z 


AB  ’ 

is  an  identity  in  the  Born  rate-of-straln  tensor 
(A.  3) 

Set 

(A.  4)  £(E)  •  E, .  £(h._)  +  K 

hAB  AB 

where  K  is  defined  for  any  P  be  formal  expansion  of  the  left-hand  side  of  (A.  4). 
If  we  decompose  the  product  in  (A.  4)  with  the  projection  operator  r,  we  have 
(A.  5)  £(E>  •  t(E,h  )•  r(£hAB)+ K+ 41 -*)(E,h  )  •  (1  -  »>UhAB>  , 

so  that  (A.  2)  and  (A.  3)  lead  to 
(A.  6) 


AB 


2  {r(E,.  )-ehoAB}«._  +  K 

nAB  AB 


♦  (1  -  s)  (£,  h  )  •  (1  -»)(£  hAB>  *  0  . 
A  B 


AB,  it  must  hold  for  Born  rigid  motions  (that  is, 


If  (A.  6)  is  to  be  an  identity  in  t 
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motions  for  which  >0).  Hence  we  obtain  the  condition 


(A.  7) 
Now, 


K  +  (1-*)(E, .  )  •  <1  -»)<£hAn>  »  0 

hAB  AB 


<!-,)(£  hAB>  ■ 


while  (l  -  »)  (E,  )  admits  the  representation 

hAB 


with 


(1-irME,  )  -  2  A^  WB)  +  B  WB 
nAB 


B  *  W.  W  ( 1  —  w>  (E, .  ) ,  AA  -  W_  (1  -  i)(E,  J-BW*. 

A  B  hAB  B  hAB 

Hence  we  have 

(1  -  x)  gE,  h  )  •  (1  -  w)  <£  hAQ)  -  2AAWA 
AB 

on  noting  that  WB  ■  0  as  a  consequence  of  WA  ^  “  1.  Combining  this 
with  the  above  result  shows  that  a  necessary  condition  for  the  existence  of  6  and 
E  is 

(A.  8)  K  +  2  AA  WA  -  0  . 

If  this  condition  is  satisfied,  (A.  1)  and  (A.  6)  imply  that* 

(A.  9) 


8  hoAB  -  »(E, .  ). 

°AB 


AB 

Conversely,  if  a  is  defined  in  terms  of  E  and  8  by  (A.  9)  and  (A.  8)  is 
satisfied,  then  (A.  4)  reduces  to  an  identity  in  <Ag*  Hence,  a  necessary  and 
sufficient  condition  for  the  existence  of  functions  8  and  E  such  that  (A.  2)  be 
an  identity  in  t  A0  Is  that  (A.  8)  and  (A.  9)  be  integrable  for  E  . 


*Since  the  rank  («A0)  $  3,  (A.  6)  Implies  8h  •  »<E,  )  +  M^, 

AB  AB 

where  « A0  *  0  .  The  condition  (A.  1)  together  with  the  fact  that  (A.  8)  is 

to  be  an  Identity  in  <AB>  and  hence  must  hold  for  «A0  of  maximum  rank,  then 

yields  »  0. 
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In  view  of  the  above  results,  we  would  like  to  establish  conditions 
under  which  (A.  9)  holds.  The  first  thing  we  must  do  is  to  bring  the  right-hand 
side  of  (A.  9)  into  a  more  tractable  form.  From  the  definition  of  kAg,  namely 


kAB  '  ,<hAB)'  tove  hAB 


+  WA  Wg.  Hence,  since  the  projectors 


B  B  B 

»A  *  &A  -  WA  W  are  idempotent  and  kAg  is  a  projection,  we  have 


AB 


'A  *B  kCD  *  W.  WB  • 


which  describes  the  projective  decomposition  of  hAg.  From  this  equation  we 

formally  obtain 

h  »  »C  eD 

nAB,kCD  *A*B- 

Hence,  if  we  replace  all  functions  of  hAg  appearing  in  L  Jby  kAg  +  WA  Wg, 
we  have 

E’k  *  E’h  hCD  ir  *  *A  E’h  *  »(C’h  1  ' 

KAB  nCD  CD’*AB  A  B  "CD  hAB 

The  system  (A.  9)  may  thus  be  written  in  the  equivalent  form 
(A.  10) 


0  h  oAB  »  L, 


AB 


The  tntegrability  conditions  for  such  systems  are  well  known**  and  all  that  must 
be  done  is  to  translate  these  conditions  into  a  form  appropriate  to  the  problem  at 
hand.  Define  the  quantities  as  follows: 


*  CD 

This  is  seen  from  the  fact  that  hAg  »  »A  »g  G^g  +  WA  Wg  gives  a 

representation  of  hAg  for  any  symmetric,  nonsingular  tensor  GAg  . 

** 

Forsyth,  A.  R. :  Theory  of  Differential  Equations.  Vol.  1,  Dover 
Publications,  Inc.,  New  York  (1959),  p.  12. 
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(A.  11) 


E ABCD(o)  >  (hflAB)„ 


(hoCD),, 


CD  ~AB 

We  then  have  the  following  result.  Ne«„  ary  and  sufflc|ent  condition*  for  the 

existence  of  functions  8  and  E  that  sail*  ^  are 

(A.  12)  eABCD  0FG  +  eCDFG  gAs  +  eFGAB  CD,  Q 

It  is  also  to  be  noted  that  the  above  condm^g  re<juce  to 
(A.  13)  gABCD  ,  o 

for  the  case  8  »  1  . 


